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ABSTRACT 

In this work we analyze a sample of AGN spectra, selected from the 6'^ Data Release of the Sloan 
Digital Sky Survey, exploiting a generalized technique of line profile analysis, designed to take into 
account the whole profiles of their broad emission lines. We find that the line profile broadening 
functions result from a complex structure, but we may be able to infer some constraints about the 
role of the geometrical factor, thus improving our ability to estimate AGN properties and their relation 
with the host galaxy. Our results suggest that flattening and inclination within the structure of the 
Broad Line Region (BLR) must be taken into account. We detect low inclinations of the BLR motion 
plane with respect to our line of sight, typically i < 20°, with a geometrical effect which generally 
decreases as the line profile becomes broader. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert — line: profiles — quasars: 
emission lines 



1. INTRODUCTION 

During the past decades, spectroscopic observations 
provided a fundamental starting point for our under- 
standing of the physical processes in Active Galactic 
Nuclei (AGN). In particular, the study of broad emis- 
sion lines, characterizing the spectra of many objects 
at different wavelengths, is a key feature t o penetrate 
their intrinsic properties (jOsterbrock I Il989l ). Unfortu- 
nately, a correct interpretation of the observations re- 
quires to collect information in various spectral ranges 
and for long monitoring times, a task which is cur- 
rently possible only in a fairly small number of cases. 
While much work has been devoted to calibrate empir- 
ical methods whi ch should be able to deal with larger 
samp les (see e. g. iKaspi et al. | [2000. 2005; Bcntz ct al. 
|2006| ). the fundamentally unknown structure in the core 
of AGN influences the physical interpretation of spec- 
tra with a geometrical factor, whose value depends on 
the structure and orientation of the sou rce (see for in- 
sta nce Vestergaard ct al. 2000; Nikolai uk et al. 1 120051 : 
iMarziani et aL II2006I : iDecarh et al. H20'08ah . 

The prominent broad emission lines, visible in the spec- 
tra of many AGN, originate close to the central power 
source, in the so called Broad Line Region (BLR). Be- 
cause of its small distance from the power source, the 
BLR is in strong interaction with the radiation field pro- 
duced by the central engine and with its gravitational 
forces. Many interesting details about the physics of pro- 
cesses that are taking place within AGN can be identi- 
fied in the signal of the BLR, but they suffer from a still 
missing complete picture of the complex kinematical and 
thermodynamical properties of the line emitting plasma. 
Since it is not yet possible to directly observe the spatial 

^ Department of Astronomy, University of Padova, 
Vicolo dell'Osservatorio, 1-35122 Padova, Italy; 

giovanni.lamura@unipd.it, stefano.ciroi@unipd.it, 
francesco.dimille@unipd.it, piero.rafanelli@unipd.it. 

^ Astronomical Observatory, Volgina 7, 11060 Belgrade, Serbia; 
lpopovic@aob.bg. ac.yu. 

^ Isaac Newton Institute of Chile, Yugoslavia Branch, 11060 
Belgrade, Serbia. 



distribution of the broad line emitting medium, although 
many important achievements were obtained in the an- 
g ular resolution of AGN cores at radio wavelengths (e. 
g. iKellerman et al~lll998l . and references therein), spec- 
troscopic data are still the most useful way to investi- 
gate physics within the BLR. The well known Reverbera- 
tion Mapping (RM) technique (sec Blandford & McKee] 
I1982D . based on multiple spectroscopic observations, pro- 
vides a reliable way to constrain the volume where 
the nuclear activity is confined and, th us, to estimate 
the mass concentration therein (e. g. IWandel | Il999t 
[Peterson & Wandel 2000; Peters on et ani2004f ). 

The most commonly accepted interpretation of AGN 
puts a Super Massive Black Hole (SMBH) in the role of 
the central engine, since matter accretion into its gravi- 
tational field provides the required power to account for 
many observational properties. Once the size of the BLR 
is known, it is possible to estimate the mass of the SMBH: 



BH 



f 



G 



(1) 



where G is the gravitational constant, / represents the 
geometrical factor which accounts for the unknown dis- 
tribution of the line emitting material, Rblr is a char- 
acteristic BLR radius, and Aw is an estimate to the ve- 
locity field, usually coming from t he width of th e emis - 
sion lines. In their recent work, iKaspi et al. I ()2005f ) 
found that a power law relationship of the form Rblr oc 
(5100 AL5100)" may adequately describe AGN, although 
the actual value of the power law exponent and the lu- 
minosity range, where the relationship hol ds in different 
AGN classes, still have to be constrained (|Kas pi et al. I 
2007tlKellv et al. l[2(M[Laori[2?M 2008). 
Bentz et al. I ()2006[ ) argued that, at least for moderate 



luminosity sources, with Lboi < 10^^ erg s^^, it is hkely 
that a ~ 0.52, not far from the predictions of simple 
photoionization calculations giving a ~ 0.5. 

Unfortunately, the problems introduced by our limited 
knowledge about the actual structure of the BLR badly 
affect the value of such estimates, rising many uncer- 
tainties which make it quite difficult to infer the physical 
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properties of AGN or to study their relation with the 
host galaxy environment. During the past years, a lot 
of work has been devoted to understand the relation- 
ship between the BLR dynamics and the correspondin, 
broad emissio n line p r ofiles ( e. g.lCap riotti et al. "f98 
1981t [F crland et al. I I l992t iPeterson & Wandcl 1999; 
Korista fc Goad lipoid but, while the former is proba- 
bly very complex, often with evidence for multiple com- 
ponents (jPopovic et al. |[2i304f) . the latter is the result of 
a combination of effects involving the gas motion pattern 
and the radiation transfer across an environment which 
is only approximately understood. 

In this paper we describe the results we obtained by 
analyzing the kinematical properties of the BLR gas in a 
sample of AGN extracted from the Sloan Digital Sky Sur- 
vey (SDSS) database, by means of a technique exploiting 
the cross-correlation method and the Gauss-Hermite pro- 
file fitting to infer the line broadening function (BF) in 
the optical domain. We show that some interesting clues 
to the geometry of the BLR can be identified in this way 
and we apply the results to estimate the physical prop- 
erties in our sample of AGN. 

The paper is organized as follows: in §2 we describe 
the analytical formalism to extract the profile broaden- 
ing functions and to calculate the corresponding Gauss- 
Hermite expansions; in §3 we present our sample and the 
reduction techniques that we adopted; §4 summarizes our 
results, with a discussion of the main limits and some in- 
dications to improve the analysis; finally our conclusions 
are given in §5. 

2. LINE PROFILE ANALYSIS AND THEORETICAL 
MODELS 

In the effort towards revealing the intrinsic properties 
of AGN cores, line profile analysis often played a ma- 
jor role. Under specific assumptions about the dynami- 
cal conditions within the BLR, such as the hypothesis of 
virial motions driven by the combined effect of the central 
engine's gravity and radiation pressure, a number of rep- 
resentative parameters, like the line widths at different 
intensity levels, or the line asymmetry factors, usually 
computed in the form of ratios among the line exten- 
sion toward the blue and red wavelengths, with respect 
to the line core position, were used in order to describe 
the profiles and to evaluate the properties of the engine. 
This kind of approach is prone to the effects of the sub- 
stantially unknown BLR geometry, with the possibility 
to introduce systematic misinterpretation of data. Fur- 
thermore, it assumes quite specific measurements to be 
a good approximation of the entire emission line profile, 
loosing some precious physical details. 

In this section we describe a generalized approach to 
the line profile fitting, already exploited in the past years 
in the field of advanced stellar kinematics, but adopted 
for gas kinematics as well (Barton et al. 2000). 

2.1. Emission line broadening from cross- correlation 

The BLR spectrum shows several emission lines 
corresponding to many permitted and some semi- 
forbidden transitions of variously ionized atomic species. 
There are indications that the distribution of the 
line emitting material is different, according to the 
ionization po tential of the con s idered emission lines 
(see e. g. iGaskell fc: Sparke! 119861 : ISulentic et al. I 



1 19951 iMarziani et al I Il996t ISnedden fc GaskeHl 
pM; Matsuok a et al. 1 120081: IMullanev fc Ward 1 120081 : 
iSlusc et al. 200l, etc.). Indeed, the interaction of 
the most energetic AGN radiation with gas probably 
produces a region where matter is highly ionized. On 
the contrary, optical shielding effects allow for the 
survival of low ionization species, in regions where 
only comparatively low energy photons may penetrate. 
Therefore, if the BLR structure is such that the shielded 
component is different from the directly exposed one, 
the properties of the emission lines will depend on their 
ionization potential. 

On the other hand, choosing to analyze a set of emis- 
sion lines belonging to a statistical distribution of matter 
and radiation interactions, it is more likely that the emis- 
sion regions are not dramatically different. In the opti- 
cal domain, the Balmer series of hydrogen is the most 
appropriate choice, because of its strength above the un- 
derlying continuum. 

Assuming that the Balmer line emission is not af- 
fected by large variations across the BLR and introduc- 
in g the cross-correla.tion formalism, originally described 
bvlTqnrv fc Davis"! ()1979D and then updated by Statler] 
()1995f ). we can approximate the observed line spectra 
as the convolution of an appropriate template of narrow 
emission lines T{x) with the BLR broadening function 
Bix): 

S{x) ~ T{x) * B{x), (2) 

where S{x) is the observed spectrum, while x represents 
a logarithmic wavelength coordinate of the form: 

x^A\n\ + B, (3) 

such that the effect of radial velocities results in linear 
shifts along x. Eq. (2) can be explicitly written as: 



S{x) ~ / T{x)B{x-x')dx' 



(4) 



and, if we compute the cross-correlation function of the 
spectrum with the template, we find: 

X{x) = S{x) (g) T{x) = / S{x)T{x + x') dx'. (5) 



Using Eq. (4), the cross-correlation function becomes: 

x,.).//r(.,i,(.-.")d.''T(. + .',d.'. (6) 

which, upon changing order of integration, is: 

X{x)~ I [ T{x)T{x + x')dx'B{x- x")dx". (7) 



Based on the definitions of cross-correlation and convolu- 
tion, Eq. (7) approximates the cross-correlation function 
of the spectrum and the template as the convolution of 
the template auto correlation function with the object's 
BF (|Statler II1995I) : 



X{x) ~ [T{x)(g)T{x)]* B{x). 



(8) 



Since T{x) is known and X{x) is drawn from observa- 
tions, as far as the template is correct, it is possible to 
recover B(x). 

Restricting our analysis to the primary cross- 
correlation peak, which carries most of the kinematical 
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TABLE 1 
List of objects in the sample 



Name 


RA 


Doc 


Mv 


z 


FWHMh/9 


Detected 






( dd ; mm ; ss ) 






(kms"-'-) 


in Radio 


RX J0801. 5+4736 


08:01:32.0 


+47:36:16 


-23.3 


0.157 


7009 ± 336 


Yes 


RXS J080358. 9+433248 


08:03:59.3 


+43:32:58 


-24.4 


0.451 


2632 ± 222 


No 


IRXS J080534.6+543132 


08:05:34.9 


+54:31:29 


-26.5 


0.406 


3174 ± 299 


No 


SDSS J081222. 99+461529.1 


08:12:23.0 


+46:15:28 


-23.6 


0.311 


2807 ± 283 


No 


2MASSi J0816522+425829 


08:16:52.3 


+42:58:30 


-23.1 


0.235 


3768 ± 186 


No 


NGC 2639 UIO 


08:42:30.5 


+49:58:03 


-23.3 


0.305 


6046 ± 210 


No 


SDSS J085632. 39+504114.0 


08:56:32.4 


+50:41:14 


-23.9 


0.234 


2819 ± 157 


No 


SDSS J085828. 69+342343.8 


08:58:28.6 


+34:23:44 


-24.3 


0.257 


3570 ± 265 


No 


SDSS J090455. 00+511444.6 


09:04:55.0 


+51:14:44 


-23.4 


0.224 


3620 ± 142 


No 


RX J0906.0+4851 


09:06:01.3 


+48:51:48 


-24.1 


0.390 


2503 ± 100 


No 


RX J0908.7+4939 


09:08:47.4 


+49:40:07 


-24.5 


0.421 


2137 ± 71 


No 


IWGA J0931. 9+5533 


09:32:00.1 


+55:33:48 


-23.1 


0.265 


5045 ± 482 


Yes 


SDSS J093653. 84+533126.8 


09:36:53.8 


+53:31:27 


-23.1 


0.228 


4255 ± 175 


No 


FIRST J094610.9+322325 


09:46:10.9 


+32:23:26 


-24.5 


0.403 


1684 ± 235 


Yes 


KUV 09484+3557 


09:51:23.8 


+35:42:48 


-23.8 


0.398 


3244 ± 289 


No 


HS 1001+4840 


10:04:13.8 


+48:26:06 


-24.0 


0.562 


2962 ± 230 


No 


PC 1014+4717 


10:17:30.9 


+47:02:25 


-23.7 


0.335 


2583 ± 103 


No 


RX J1030.4+5516 


10:30:24.9 


+55:16:21 


-24.8 


0.435 


2246 ± 157 


Yes 


FBQS J103359.4+355509 


10:33:59.5 


+35:55:09 


-23.3 


0.169 


4543 ± 237 


Yes 


SBS 1047+557B 


10:50:55.1 


+55:27:23 


-24.1 


0.333 


2253 ± 97 


No 


RX J1054.7+4831 


10:54:44.7 


+48:31:39 


-24.7 


0.286 


4428 ± 274 


Yes 


FBQS J105648.1+370450 


10:56:48.2 


+37:04:51 


-24.1 


0.387 


4538 ± 272 


Yes 


FBQS J110704.5+320630 


11:07:04.5 


+32:06:30 


-23.3 


0.243 


8447 ± 558 


Yes 


FBQS J112956.5+364919 


11:29:56.5 


+36:49:19 


-24.4 


0.399 


3564 ± 452 


Yes 


FBQS J115117.7+382221 


11:51:17.8 


+38:22:21 


-25.4 


0.336 


3696 ± 242 


Yes 


RX J1200.4+3334 


12:00:28.7 


+33:34:43 


-23.5 


0.284 


2864 ± 173 


No 


RX J1203.8+3711 


12:03:54.8 


+37:11:37 


-23.4 


0.401 


3725 ± 1020 


Yes 


IRXS J121759. 9+303306 


12:17:59.3 


+30:33:03 


-23.8 


0.363 


3435 ± 258 


No 


RX J1218.3+3850 


12:18:22.6 


+38:50:43 


-24.4 


0.194 


4564 ± 483 


No 


FBQS J122035.1+385317 


12:20:35.1 


+38:53:17 


-25.3 


0.376 


1899 ± 164 


Yes 


FBQS J122424.2+401510 


12:24:24.2 


+40:15:10 


-24.3 


0.415 


4721 ± 397 


Yes 


FBQS J122624.2+324429 


12:26:24.2 


+32:44:29 


-23.9 


0.242 


4535 ± 533 


Yes 


FBQS J125602.0+385230 


12:56:02.0 


+38:52:30 


-24.1 


0.419 


1978 ± 147 


Yes 


FBQS J132515.0+330556 


13:25:15.0 


+33:05:56 


-23.5 


0.356 


1921 ± 102 


Yes 


SDSS J144050.77+520446.0 


14:40:50.7 


+52:04:46 


-23.4 


0.320 


2483 ± 134 


No 


RX J1452.4+4522 


14:52:24.7 


+45:22:24 


-24.6 


0.469 


6594 ± 156 


Yes 


FBQS J145958.4+333701 


14:59:58.4 


+33:37:01 


-25.5 


0.644 


4389 ± 209 


Yes 


FIRST J154348. 6+401324 


15:43:48.7 


+40:13:25 


-23.4 


0.318 


5227 ± 391 


Yes 


SDSS J154833. 03+442226.0 


15:48:33.1 


+44:22:26 


-23.5 


0.322 


3893 ± 343 


Yes 


FBQS J155147.4+330007 


15:51:47.4 


+33:00:08 


-24.5 


0.422 


2744 ± 221 


Yes 



Note. — Together with names and equatorial coordinates (J2000.0), the table lists the absolute magnitude in the V band, the 
cosmological redshift, FWHMh^, and a binary flag indicating whether the source is detected in radio or not. 



information and it is weakly affected by template mis- 
match, Eq. (7) can be written in its discrete form, with 
the simplified notation Fi ~ F{xi): 



N 



N 



i+j 



(9) 



1=0 \i=o 



Provided that all the functions arc null when they are 
computed outside the range < i < iV, Eq. (9) defines a 
system of iV + 1 linear equations in the TV -I- 1 variables 
Bk~^i {k > i). A standard -minimization routine can 
be therefore used to infer the BF of the Balmer lines. 

2.2. Analytical expressions for the broadening functions 

As previously mentioned, the BLR broadening func- 
tions are influenced by the effects of complex kinemat- 
ics within the source and of radiation transfer from the 
source to the observer. For this reason it is hardly con- 
ceivable that a simple analytic expression might be used 
to fit the resulting profiles. In the case of a geometrically 
complex distribution of motions in the line emitting re- 
gion, multiple Gaussian functions provide reasonable fits 



to the observed profiles. Two Gaussian c ontributions can 
usually fit the broad co mponent of H/3 (jPopovic et al. I 
12004 IChen et al. l[2008l ). but other contributions, up to 
five more Gaussians, might be needed to account for the 
narrow emission lines of H/3 and [O III]. Furthermore, 
the presence of ordered kinematical components modi- 
fies the shape of the BF, raising non-Gaussian features 
in the profiles. 

A good way to estimate the importance of non- 
Gaussian components is to parameterize the ob- 
served BF by means of a Gauss-Hermite orthonor- 
mal expansion, similar l y to what is described in 
IVan Per Marel k. Franx I (|1993f) for the case of stellar 
kinematics in elliptical galaxies. Following their method, 
if we call a{v) the normal Gaussian function: 



a{v) 



cxp 



2al 



(10) 



where is the line of sight velocity dispersion, the emis- 
sion line BF can be expressed as a function of v : 



B{v) = B^aiv - Vsys) 



N 



l + ^h,H,{v-Vs, 



, (11) 
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in which we caU Bq the BF normahzation factor, Vgys the 
systemic radial velocity offset between the BF and the 
chosen reference frame, Hi(v—Vsys) the i**^ order Hermite 
polynomial, and hi the corresponding coefficient. A wide 
descriptio n of the properties of the Hermit e polynomials 
IS given m iVan Per Marel fc Fraiixl ()1993f) . It is demon- 
strated that odd order functions account for asymmet- 
ric distortions of the Gaussian profile, while even order 
functions have a symmetric effect. Truncating Eq. (11) 
to iV = 4, the Hermite polynomials are expressed by: 



i/3(2/) = ^(2V2y3_3V22/) 



(4y4_i2y2 + 3). 



(12a) 



(126) 



Therefore, it is possible to estimate the role of non- 
Gaussian kinematical components, using the whole BF 
profile, simply by fitting the observed shape with a trun- 
cated Gauss-Hermite series and measuring the appropri- 
ate values of and /14. 

3. SAMPLE SELECTION AND DATA REDUCTION 

3.1. AGN spectra from SDSS 

To perform our investigation, we needed a homoge- 
neous sample of AGN optical spectra, featuring promi- 
nent broad Balmer emission lines. The spectroscopic 
database of the 6*^ data release of SDSS (DR6) provides 
a huge number of objects, whose spectra are collected and 

S irocessed by a fairly well establi shed standard pipeline 
Adelman-McCarthv et al. Il2008f) . In order to collect a 
set of good signal spectra, we chose to select the sample in 
the Veron-Cetty cata logue of Quasars and Active Galac- 
tic Nuclei (12*'* Ed.) ()V^ron-Cettv k Veron |[2006h . with 
the following requirements: 

• the object redshift had not to exceed the limit z k, 
0.8, because objects at larger redshift have their H/3 
emission line beyond the SDSS spectral coverage; 

• only bright sources, with My < — 23, were consid- 
ered; 

• each object had a spectral classification suggesting 
the presence of broad emission line profiles; 

• the sources are located within the field covered by 
the SDSS DR6 spectroscopic observations. 

On the resulting candidate list, we applied further con- 
straints, that restricted the sample to the most appropri- 
ate sources. In particular we chose objects whose spectra 
had at least three clearly detectable broad Balmer lines 
and they were not affected by instrumental disturbances 
or by strong foreground and background contamination. 
We ended up with a sample of 40 objects that is described 
in Table 1. 

3.2. Preliminary spectral reduction 

A major advantage in the SDSS database is that it 
provides spectra with preliminary reduction and cali- 
bration, thus simplifying the task of spectral analysis. 
Therefore, before proceeding with our measurements, 
we simply had to remove from the spectra those con- 
tributions which come from outside the AGN. We ap- 
plied a correction for Galactic Extinction, estimated by 
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Fig. 1. — Example of spectral decomposition for SDSS 
J085632. 39-1-504114.0. The observed spectrum, represented by the 
thick continuous line, is compared with the AGN (thin continuous 
line) and galactic (long dashed) components. The dotted line in 
the bottom part of the plot is the fit residual. 

means of a selective extinction f unction, in the form pro- 
posed by ICardelli et al. I (|1989l ) with the absorption co- 
efficients evalu ated on the basis of the extinction map 
traced by Schlcgcl et al. I (|1998) and available at the 
NASA Extragalactic Extinction Calculator. We, then, 
removed the cosmological redshift, bringing the spec- 
tra to the rest frame of the narrow [O III] emission 
lines. Finally, we estimated the host galaxy contami- 
nation by applying a spectral decomposition technique, 
based on the Karhune n-Loeve Transforms described by 
IConnollv et al. I (I1995D a,nd im plemented on SDSS data 
by IVanden Berk et al. I (I2006D . According to the same 
method exploited in iLa Mura et al."l (j2007l ). we consider 
the observed spectra as the linear combinations of princi- 
pal components, called eigenspectra, originated indepen- 
dently by the AGN and its host: 



^(A)=^fe-Q.(A)]-l-^[g,-G,(A)], 



(13) 



with S'(A) being the total spectrum, Qi{X) the z*'' AGN 
component, weighted by its coefficient qi, and Gj{X) 
the j'^ host galaxy eigenspectrum associated to the 
corresponding coefficient gj. Using the galact ic and 
AGN eigenspectra provided bv lYip et al. I (|2004al lbl) . we 
evaluated the coefficients (g^, gj) by means of a x^- 
minimization routine involving the first five galaxy eigen- 
spectra and six AGN components. This procedure allows 
to carry out the separation outlined in Eq. (13), as it is 
illustrated in Fig. 1, and to subtract an estimate of the 
host galaxy starlight contaminating the SDSS spectra. 

3.3. Extracting the BLR component 

Once the AGN spectra have been corrected to account 
for most of the external effects, the task to identify the 
BLR contribution alone needs the removal of more com- 
ponents, including the underlying continuum of the AGN 
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Fig. 2. — Multiple Gaussian decompositions of the profiles of H/3 (upper left panel), Ho (upper right), H7, and H(5 (lower panel). Here 
we use a thick continuous line to plot the continuum subtracted spectrum of PC 1014+4717, a thin continuous line for the estimated NLR 
contributions, a long dashed line for the BLR components, and a dotted line for Fe II. The thick continuous line in the bottom part of each 



panel shows the fit residuals, while fluxes and wavelengths are expressed in units of 10 ^' erg cm s 



and A, respectively. 



central source and the narrow lines coming from the Nar- 
row Line Region (NLR). Measuring the properties of the 
broad lines, moreover, is often difficult because of the 
multiple spectral features that are blended together. In 
the case of the Balmer series, narrow lines from [O III], 
[N II], [S II], together with the narrow Balmer emissions, 
have to be taken into account, while blends with broad 
lines from He II and the multiplets of Fe II can heavily 
affect the observed profiles. 

To subtract the AGN continuum, we fit the spectra 
with spline functions of order ranging from 2 to 5 
in wavelength ranges which are usually not affected by 
prominent lines. The subtraction of narrow lines and 
the blends with He are dealt with by means of multiple 
Gaussian profile decompositions, where we use the [O III] 
emission line at 5007 A as a template for the other nar- 
row features and we fix the [N II] AA 6548, 6584 emission 
line ratio to be 1:3. Some special care, instead, is needed 
in the case of the Fe II multiplets, whose properties have 
not yet been understoo d in detaiL Follo wing the method 
suggested bvlVeron-Cct tv et al. I ()2004D . it is possible to 
remove the Fe II contribution from spectra by scaling and 
broadening an appropriate template, estimated from the 
spectrum of an AGN featuring prominent Fe II emission 
lines. In our work, we used the Fe II template spectrum 



coming from / Zwicky 1 (jBotte et al. II2004D . splitting it 
into two parts, which we scaled in order to achieve a bet- 
ter coincidence with our data. The process is summarized 
in Fig. 2. Most of the steps so far described were carried 
out with tasks provided in the IRAF software package. 

Once the removal of contaminating contributions has 
been performed, we are able to apply the Gauss-Hermite 
formalism to the profiles of the broad lines so far isolated 
in the spectra. We perform a first analysis directly on 
the profiles of the H/3 emission line. The results of this 
study can be subsequently compared with the shape of 
the Balmer line broadening functions, which we are now 
able to infer from the cleaned BLR spectra. 

3.4. The Balmer Line Broadening Functions 

As a result of the previous steps, we now have got 
a set of BLR Balmer line spectra. Our task is then 
to recover their BF, by means of the cross-correlation 
technique outlined in §2.1. To calculate the cross- 
correlation functions, we build a template of Balmer 
emission lines, following the me dian line intensity ratios 
found bv lLa Mura et al. I (|2007[) in a sample of 90 SDSS 
spectra of broad line emitting AGN. Our template as- 
sumes that the SDSS instrumental profile is a Gaussian 
function with FWHM = 167 km s^^. At the spectral res- 
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Fig. 3. — Gaijss-Hermite expansion coefficients of the line profile BF compajred with the same coefiicients extracted from the profile of 
H/3. The straight lines shown in both panels axe the best fit functions described in Eq. (16) and Eq. (17). 



olution of Sloan data, the logarithmic sampling of the 

wavelength coordinate can be performed with discrete 
bins corresponding to 69 kms^^ each. Here we use the 
IRAF task f xcor to compute the template autocorrela- 
tion function: 

A{x) = T{x) T{x) (14) 
and the cross-correlation functions of the BLR spectra 
with the template, following the definition of Eq. (5). 
Again a x^-minimization algorithm can be exploited to 
infer the BF in its discrete form. Applying the least 
squares formalism to the equation system (9), it follows 
that the BF of each spectrum must satisfy the relations: 

N / N \ 

Y.bAJ2 = J2 Ai-kXi. (15) 

i=0 \j=0 J i=k 

In principle, it is possible to extract an accurate solu- 
tion for the BF by solving the equation system (15) with 
< fe < A'^. In practice the task is not simple, because 
it involves the inversion of a coefficient matrix as large 
as [(A^ + 1) X (Af -I- 1)], with N increasing with the line 
profile widths up to A" ~ 400. However, the complete 
solution of such a system is not the real purpose of this 
work, since we arc not seeking the detailed shape of the 
BF, but we are rather looking for the importance of non- 
Gaussian components. Therefore, we chose to solve the 
system at lower resolution, interpolating the BF every 8 
bins with an analytical profile, which we assumed to be 
a Gauss-Hcrmite expansion. 

We compared the properties inferred for the BF of 
our spectra with the results obtained by applying the 
Gauss-Hermite profile fitting directly to the H/3 emission 
line. We found that the expansion coefficients in the two 
cases are highly correlated, supporting a tight relation- 
ship among the H/3 emission line profile and the BF of 
the Balmer series. The details of this comparison are 
given in Fig. 3, where we plot the values of the expan- 
sion coefficients obtained in both ways. As a result, we 
get: 

/if ^ = (0.773 ± 0.073)/i3 '^^ + (0.721 ± 0.554) • 10-^ (16) 
with a correlation coefficient R — 0.865 and a null hy- 
pothesis Pq < 10~^ and 

/if ^ = (1.051 ± 0.052)/i4 + (0.154 ± 0.049) • 10"^, (17) 



with R = 0.955 and Pq < 10"^ 

Here we would like to point out that the matrix in- 
version must be computed only once, because it involves 
coefficients exclusively drawn from the template autocor- 
relation function. The cross-correlation functions of the 
spectra, instead, only affect the known terms of Eq. (15). 
Hence the advantage of this technique. 

3.5. Spectral property rn,easurem,ents 

After the calculation of the BF of our sample, we per- 
formed more measurements of spectral properties in the 
data, estimating, in particular, the FWHM of the H/3 
emission line and the AGN continuum radiation lumi- 
nosity at 5100 A. These are needed to infer some of the 
source physical properties, such as the central black hole 
mass, its accretion rate and the size of the BLR. 

In the case of the line profile measurements, we looked 
at the H/3 emission line in the BLR spectra, which we pre- 
viously isolated for cross-correlation with the template. 
To account for the uncertain continuum and narrow line 
corrections, we performed five different measurements of 
the line half width at half the maximum both on the blue 
and red wing of the line, varying our guess to the con- 
tinuum and line peak intensity. We combined these es- 
timates to calculate the FWHM, then we averaged them 
together, and we computed the la dispersion. As a fur- 
ther step, we fit the broad H/3 profile with two Gaussian 
functions and we applied similar measurements to the 
identified components. 

The optical continuum luminosities at 5100 A, instead, 
were evaluated in the AGN spectra that we previously 
corrected for Galactic Extinction and host contamina- 
tion. Here, the main source of error arises from the noise 
fluctuations around the actual signal intensity. For this 
reason, we assumed the specific continuum luminosity 
at 5100 A to be representcid by the average luminosity, 
evaluated in the range running from 5075 A to 5125 A, 
and the associated error to be given by its standard de- 
viation. Therefore, we measured the continuum fluxes 
of our spectra and we computed the related specific lu- 
minosities using the object redshift as a distance indica- 
tor, in the framework of a cosmological model defined by 
Ho = 75kms~^ Mpc~^, flmatter = 0.3, and Oa = 0.7. A 
guess to the bolometric luminosity, which is needed to in- 
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Fig. 4. — Expected FWHM in the BP of emission lines origi- 
nated in a disk structure. The empty triangles connected by the 
continuous line are the predicted values for reference disk models 
with i = 10° in the gravitational field of black holes with M = 1, 
2, 4, 6, 8, 10, 12, and 14 ■1O''M0 (respectively models from 1 to 
8); the circles with the dashed line show the predicted behavior for 
the reference disk model in the gravitational field of the black hole 
with M = IO^'Mq, when seen under inclinations of i = 10°, 15°, 
20°, 25°, 30°, 35°, 40°, and 45° (again from model 1 to 8). 

fer the accretion rate onto the central black hole, can be 
made from the specific luminosity measured in the spec- 
tra (e. g. lElvis et al. I '1994: b ut see als o Collin et al. I 



l2002i ICollin fc Kawaguchi 2004 and references therein) . 
Here we assume our objects to have bolometric luminosi- 
ties approximately given by: 



^bol 



10- (5100A-L5ioo), 



(18) 



where L5100 represents our estimate to the continuum 
specific luminosity at 5100 A, measured in ergs~^ A^^, 
and Lboi is the bolometric luminosity, given in ergs~^. 
We note, however, that this assumption is prone to the 
effect of the large dispersion in the SED of AGN and 
it may easily introduce an uncertainty of a factor ~ 2, 
which propagates in the estimated black hole masses and 
accretion rates. 

4. RESULTS AND DISCUSSION 

A clear determination of the black hole mass is not pos- 
sible unless we are able to discriminate the role played 
by / in Eq. (1). In many circumstances, the black 
hole mass problem is dealt with through the assump- 
tion of a particular geometry, such as a random distribu- 
tion of virial motions (Nctzcr 1990; Wandcl et al. 1999; 
iPeterson fc Wandel l 1200 0), or a flattened rotating sys- 
tem with an inclinatio n mostly inferred by rn eans of sta- 
tist icalconsiderationsj[Dec^^ 

thors (IVes tergaard et a l. I l2000t iNikolai uk et al. 1 120051 : 
lPeterson"e t al. 2004; S ulentic et al. 200(|, etc.) pointed 
out that a considerable flattening is likely to be an in- 
trinsic property of the BLR structure and even that a 
simple assumption about its inclination may not remove 
the problem. The very nature of the broad line emitting 
entities has been inves tigate d extensively (see, for exam- 
ple, lArav et al. Ill997l 119981: iLaor Il2006if) . leading to the 



conclusion that the broad line proflles either result from 
the combination of a large number of emitters, in the or- 
der of ^ 10^, or it is produced by motions of a smooth 
medium. In both cases, it is noticed that a random mo- 
tion pattern could not be dynamically stable. 

With the exception of those objects whose line pro- 
flles clearly show double peaks, a strong clue towards an 
highly inclined rotating system, the BLR inclination is 
still an open question of crucial importance for determi- 
nation of black hole mass and accretion rate. Indeed, 
the assumption of a universal geometrical factor usually 
leads to the detection of high accretion rates in Nar- 
row L ine Seyfert 1 galaxies (NLSl) (e. g. iBoUer et al. I 
Il996f ). while adjusting the geometrical factor, according 
to statistics, affects the black hole masses largely reduc- 
ing most of the differences. Both these paths might be 
sources of systematic misunderstandings, therefore a di- 
rect and independent measurement of the BLR inclina- 
tion, or, alternatively, of the black hole accretion rate, 
would be highly desirable in order to discriminate among 
the a ctual dynamical pr operties and the effect of inclina- 
tion (Kell v et al.1l2008D . 

A si milar test has beeri perfor med, for a restricted sam- 
ple, bv lHicks &: Malkan I |2008h in near infrared spectro- 
scopic observations, leading to the conclusion that the 
observed gas kinematics is consistent with RM based re- 
sults for nearly face-on disk structures. 

4.1. Inclination and line profile broadening 

Since the BLR structure cannot be represented by a 
random motion pattern, the shape of the broad emission 
lines exhibits large deviations from the Gaussian proflle. 
In the extreme case of the marked do uble peaks in the 
spectr al lines of Arp 102B, studied bv lChen fc Halpern I 
()1989t ). the BLR structure is well explained in terms of 
the combination of a quasi spherical component with a 
rotating disk, probably the external accretion disk, seen 
at an inclination of i = 32°. Popovic ct al. (2004) ap- 
plied the same model to other single peaked line emit- 
ting sources and they found that mildly inclined disks 
could be responsible for the observed line proflles as 
well, although uncertainties on the model free param- 
eter s might affect the values of the inferred inclinations, 
as iCollin et al. I (|2006f ) pointed out. 

If the BLR has a flattened component which is seen at 
low inclination, its emission lines clearly do not exhibit 
double peaks, but the geometrical structure still modifles 
the dynamical interpretation of data. We illustrate this 
concept in Fig. 4, where we plot the expected FWHM in 
the broadening function of disks surrounding black holes 
of increasing mass and we compare it to the situation of 
a black hole of flxed mass, but with the disk seen under 
different inclinations. It is clear that, within the range of 
our calculations, there is a mass - inclination degeneracy 
on the resulting FWHM. 

Exploiting the m odel develo ped inlChen et al. I (I1989D . 
IChen fc Halpern I (|1989D . and lPopovic et al. I (|2004D . we 
computed a range of expected non-Gaussian profiles in 
the case of a two component BLR structure seen at dif- 
ferent inclinations, with a flattened rotating disk and 
a surrounding distribution of gas, giving rise to a bell- 
shaped contribution. In its original purpose, this model 
was conceived to flt the properties of an accretion disk, 
introducing some free parameters for the disk radii, in- 
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Fig. 5. — Comparison of the reference model predictions con- 
cerning the broadening function kurtosis (continuous lines) with 
four variants: in the upper panel we plot models with stronger 
(long dashed line) and weaker (short dashed line) disk emission 
with respect to the bell-shaped component; in the bottom panel 
we show the differences obtained by setting crg^u = 0.07c (short 
dashed line) and crg^n = 0.09c (long dashed line). 

trinsic velocity dispersion, and line emission. Adjusting 
these parameters, it would be possible to fit the broad 
emission line profiles of possibly all the spectra of our 
sample, but reasonable fits can be obtained in several 
ways, without tightly constraining the physical proper- 
ties of the BLR. Here we try to predict the effect of a 
flattened BLR component on the observed line profiles, 
therefore we fix some of these parameters on the basis 
of the results collected by [La Mura et al. (2007). Our 
reference model assumes = 1834i?5 for inner radius, 
Rblr = 18340i?s for outer radius, aoisk = 0.003c as the 
intrinsic velocity dispersion in the disk, and a = —2.0 for 
the radial emission power law, where Rs and c represent 
the Schwarzshild radius and the speed of light. More- 
over, this model includes a bell-shaped gas distribution 
having a velocity dispersion of aBeii = 0.008c. The ref- 
erence model carries out the best match to the observed 
line profiles with the assumption of various disk inclina- 
tions. 

In Fig. 5, we compare the reference model with some 
variants, obtained with slightly different parameters. We 
note that all the models predict a strong dependence of 
the line profile kurtosis (the coefficient /14 in the Gauss- 
Hcrmite expansion) on the disk inclination, in the range 
of small values of i. The reason is quite simple, because a 
nearly face-on disk enhances the low radial velocity peak 
of the BF, increasing the kurtosis of the profile, while 
an edge-on disk is more likely to affect the high velocity 
wings. However, differences in the relative normaliza- 
tion of the bell-shaped component, with respect to the 
disk, or in its intrinsic velocity dispersion, may also af- 
fect the inferred kurtosis. The assumed strength of the 
bell-shaped component has a large effect on the predicted 
kurtosis for « < 10°, while its velocity dispersion has a 
weaker influence in the range 10° <i< 20°. Since large 
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FWHMh^ (10^ km s-') 

Fig. 6. — BF kurtosis plotted as a function of FWHMh/j. We 
observe a large evolution of the line profile kurtosis as a function 
of broadening, suggesting the possibility of inclination effects in 
nearly face-on flattened structures. Such effects become weaker as 
the profile width increases. 

changes in the model parameters quickly result in pre- 
dictions that do not match the observed line profiles, we 
assume, in our calculations, a confidential uncertainty of 
Ai = ±2°, for i < 10°, and of Ai = ±5°, for i > 10°, 
where the dependence of kurtosis on inclination becomes 
shallower. At i > 20° the kurtosis is no longer a useful 
indicator of inclination. 

As we show in Fig. 6, however, where we plot the mea- 
sured values of /14 as a function of FWHMh/3 , there is a 
remarkable evolution of the line profile kurtosis, which 
decreases for increasing line profile width. Such an effect 
is a clear indication that a considerable variation of the 
geometrical factor / might be present and it should be 
taken into account in order to estimate the actual proper- 
ties of the SMBH located in the centre. Using the model 
predictions, we can exploit the broad line kurtosis to es- 
timate the inclination of the fiattened BLR component 
and to apply a correction to our dynamical interpreta- 
tion of the observed line profiles. It should be noted, 
however, that, although the kurtosis is estimated from 
the whole profile, it reduces the available information to 
a single parameter. It is, therefore, very important that 
the model provides a good fit of the observed line profiles. 

4.2. Mass and accretion rate estimates 

It can be shown that completely neglecting the role 
played by the BLR geometrical factor may lead to incor- 
rect black hole mass estimates, with uncertainties that, 
in the worst cases, could span over two orders of mag- 
nitude. This problem is particularly important in the 
case of NLSl galaxies, whose nature has been carefully 
investigated, to find out whether they are characterized 
by flattened rotating struct ures seen at low inclination 
(jOsterbrock fc Pogge I Il985( ) , or they are actually low 
mass black holes accreting at very high rates, sometimes 
well beyond the Eddington limit, as it is discussed, for 
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Fig. 7. — BLR inclination inferred by comparison among the BF 
kurtosis, predicted by the reference model for various inclinations, 
and the corresponding distribution observed in our data. Accord- 
ing to the model predictions, the BF kurtosis is very sensitive to 
flattening and inclination in nearly face-on structures, while it be- 
comes a weaker indicator for larger inclinations. 



instance, in lBoller et al. I ()1996l ) or in iKomossa I (|2008f ). 
Moreover, the role played by non gravitational forces, 
especially in the case of high radiative efficiency, may 
also influ ence the kinema t ical p roperties of gas, as sug- 
gested bv iMarconi et al. 1 ()2008l ). affecti ng the reliability 
of the virial assumption. In their work, iLa Mura et al. I 
(|2007!) found that, although NLSl had quite high ac- 
cretion rates, they were not exceptional with respect to 
other AG N in the sample , a resul t echoed by the consid- 
erations of iDecarli et al. I (l2008aD. 

On the other hand, while emmer et al. I ((20061) ar- 
gue that X-ray observations may provide a direct clue to 
the black hole accretion rate and thus remove the degen- 
eracy i ntroduced by the FWH MHa dependent ma s s esti - 
mates, iDecarh et al~l ()2008bD and iLabita et al~l (|2006f ) 
use the black hole correlation s with the ho st properties, 
identified bv iFerrarese et al. I (|2000l . l2006f) . to calibrate 
the geometrical factor. Both methods suggest that some 
care should be taken in using only the profile of H/? to 
infer the physical properties of AGN. 

With the information coming from the line profile dis- 
tortions, we compute our estimates of the black hole mass 
and accretion rate by introducing an equivalent velocity 
field, defined by: 



%/3 



FWHMBen(H/?) 



FWHMc.,.fc(H/3) 
4sini 



(19) 

Assuming that the line profile broadening results from 
both planar and non-planar motions (J ^abita et al. 200(|; 
iMcLure et all 120021 : [Jarvis fc McLure 1 12006 , etc.), v^g 
combines the velocity estimates obtained from the H/3 
emission line profile by fitting two Gaussian functions, 
which are subsequently compared with the reference 
model, providing a distinction among the bell-shaped 
and the flattened contributions. The corresponding ge- 



ometrical factors are ass umed to be, re spectively, the 
classical interpretation of 'Netzer ' (f 990) and that of a 
rotating disk, confined in a smaller region with respect 
to the other component. The inclination of the disk is 
estimated by comparison of the BF kurtosis with that of 
the reference model, as shown in Fig. 7, and its charac- 
teristic radius is assumed to be approximately four times 
smaller than the typical size of the bell-shaped compo- 
nent. To calculate the black hole mass, we introduce 
Veq in Eq. (f), bringing the geometrical factor into the 
modified velocity field, and we estimate the correspond- 
ing Eddington ratio from the bolometric luminosity in 
Eq. (f8). 

The results of our measurements and calculations are 
summarized in Table 2, together with the computed un- 
certainty ranges. Our method breaks the strong depen- 
dence of the black hole accretion rate and mass estimates 
on the width of H/3 since it exploits more indications, 
coming from the broadening functions of the observed 
Balmer lines. Moreover, the typical inclinations inferred 
for the BL R of our spectra are co nsistent with those es- 
timated bv iPopovic et al. 1 (|2008f ). suggesting that this 
situation is quite common in single peaked broad line 
emitters. Table 2, however, does not include the errors 
which could be introduced by our assumptions, concern- 
ing the source luminosity and the BLR structure. Such 
uncertainties may be as large as a factor of 2 or 3, as 
seen in the AGN SED distribution, or by adopting dif- 
ferent structural models for the BLR. We shall further 
discuss the role of these uncertainties with the help of 
some consistency checks. 

4.3. Discussion 

While our estimates of bolometric luminosity and, con- 
sequently, black hole mass and accretion rate are es- 
sentially scaled by our measurement of the optical con- 
tinuum luminosity, the dynamical interpretation of the 
line profiles still suffers from undeniable shortcomings. 
Adopting a two component model to explain the line 
profile broadening complicates the relationship among 
FWHMh/3 and the black hole mass, introducing a ge- 
ometrical factor which depends on the inclination of 
the flattened component and on its relative importance 
with respect to the BLR as a whole. Because disks are 
the most viable solution to support accretion flows in 
presence of angular momentum, numerous authors sug- 
gested that the broad line gas could origin ate in the 
disks themselves (e. g. Shields f977; Shlos man et al. I 
Il985t lEmmering et al.1 ri992f ) . Models based on accre- 
tion disks only, however, have great difficulti es in ac- 
count ing for AGN observational properties (see iKinnev I 
11994 for example). The assumption of a two compo- 
nent model improves our ability to understand the ob- 
served line profiles, but it still fails in placing strong con- 
straints on the structure of the BLR, since the origin of 
the bell-shaped component is not clear. Indeed, there 
are rn od els, such as tho s e of | Collin-So uffrin fc Dumont J 

an£| 



1990D . iJackson et aTl ()1991[ ). or Murrav & Chi 



19971 ). which achieve a good match with observations on 



more physical grounds, either exploiting very large disk 
radii, or computing the effect of radiation transfer across 
radial gas flows close to the disk. Clearly, the choice of 
different models affects the interpretation of AGN dy- 
namical properties and this is a major concern in the 
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TABLE 2 

Estimated properties of the sources 



Name'' 




hs 


hi 


i 


Mbh 






(10'"'ergs~i) 


(10-7) 


(10-7) 


(dog) 


(IO^Mq) 





RX J0801. 5+4736* 


20 + 1 


3.1 


0.2 


18 


14 + 5 


0.011 + 0.004 


RXS J080358. 9+433248 


36 + 2 


8.9 


7.3 


12 


5 + 2 


0.054 + 0.027 


IRXS J080534. 6+543132 


83 + 2 


-71.6 


2.8 


14 


8 + 3 


0.083 + 0.032 


SDSS J081222. 99+461529.1 


37 + 1 


68.6 


8.6 


11 


10 + 4 


0.028 + 0.013 


2MASSi J0816522+425829 


25 + 1 


-0.3 


4.1 


13 


9 + 4 


0.021 + 0.010 


NGC 2639 UIO 


17 + 1 


46.8 


1.7 


15 


4 + 1 


0.031 + 0.012 


SDSS J085632. 39+504114.0 


39 + 1 


39.3 


11.2 


10 


5 + 3 


0.057 + 0.036 


SDSS J085828. 69+342343. 8 


37 + 1 


-15.4 


3.2 


14 


7 + 3 


0.038 + 0.016 


SDSS J090455. 00+511444.6 


19 + 1 


45.0 


1.9 


15 


7 + 3 


0.021 + 0.009 


RX J0906. 0+4851 


35 + 1 


197.8 


26.1 


6 


13 + 7 


0.020 + 0.011 


RX J0908. 7+4939 


87 + 2 


83.1 


25.3 


6 


14 + 6 


0.048 + 0.023 


IWGA J0931. 9+5533* 


26 + 1 


-4.7 


1.1 


16 


32 + 9 


0.006 + 0.002 


SDSS J093653. 84+533126. 8 


25 + 1 


-54.9 


0.9 


16 


36 ± 12 


0.005 + 0.002 


FIRST J094610. 9+322325* 


31 + 2 


-161.1 


21.2 


7 


5 + 2 


0.048 + 0.021 


KUV 09484+3557 


22 + 1 


49.4 


4.1 


13 


5 + 2 


0.037 + 0.017 


HS 1001+4840 


24 + 3 


58.8 


2.6 


14 


12 + 4 


0.016 + 0.008 


PC 1014+4717 


25 + 1 


85.5 


12.1 


10 


5 + 3 


0.038 + 0.023 


RX J1030. 4+5516* 


40 + 4 


144.8 


23.4 


6 


12 + 6 


0.026 + 0.016 


FBQS J103359. 4+355509* 


13 + 1 


26.7 


3.6 


14 


3 + 1 


0.030 + 0.010 


SBS 1047+557B 


43 + 2 


44.6 


18.0 


8 


8 + 3 


0.039 + 0.016 


RX J1054.7+4831* 


49 + 2 


31.1 


0.9 


16 


12 + 5 


0.030 + 0.013 


FBQS J105648. 1+370450* 


20 + 2 


19.8 


1.5 


15 


5 + 2 


0.028 + 0.014 


FBQS J110704. 5+320630* 


18 + 1 


-35.5 


1.6 


15 


24 + 11 


0.006 + 0.003 


FBQS J112956.5+364919* 


10 + 1 


67.7 


4.0 


13 


1.8 + 0.7 


0.041 + 0.020 


FBQS J115117.7+382221* 


188 ± 4 


-17.8 


5.2 


13 


22 + 8 


0.065 + 0.026 


RX J1200.4+3334 


9 ± 1 


10.8 


7.9 


11 


3 + 1 


0.023 + 0.014 


RX J1203.8+3711* 


27 + 2 


79.2 


4.2 


13 


3 + 2 


0.068 + 0.041 


IRXS J121759. 9+303306 


17 + 1 


33.4 


2.3 


14 


5 + 2 


0.025 + 0.013 


RX J1218.3+3850 


15 + 1 


-17.4 


1.9 


15 


6 + 3 


0.018 + 0.009 


FBQS J122035.1+385317* 


93 + 2 


-84.3 


23.8 


6 


9 + 4 


0.079 + 0.040 


FBQS J122424.2+401510* 


35 + 2 


15.0 


1.8 


15 


8 + 2 


0.035 + 0.013 


FBQS J122624.2+324429* 


7 ± 1 


12.7 


1.8 


15 


3 + 1 


0.016 + 0.007 


FBQS J125602.0+385230 


31 + 3 


-78.5 


22.3 


7 


4 + 2 


0.054 + 0.032 


FBQS J132515.0+330556 


27 + 1 


-29.5 


18.8 


8 


5 + 2 


0.039 + 0.017 


SDSS J144050.77+520446.0 


20 + 1 


67.7 


11.2 


10 


6 + 3 


0.028 + 0.016 


RX J1452.4+4522* 


88 + 4 


4.8 


-0.2 


19 


6 + 2 


0.105 + 0.044 


FBQS J145958.4+333701* 


89 ± 20 


41.7 


1.9 


15 


8 + 4 


0.081 + 0.056 


FIRST J154348. 6+401324* 


29 + 2 


-53.0 


3.9 


13 


5 + 2 


0.042 + 0.016 


SDSS J154833. 03+442226.0* 


23 + 1 


-6.6 


2.6 


14 


8 + 3 


0.021 + 0.009 


FBQS J155147.4+330007* 


32 + 3 


12.0 


6.4 


12 


3 + 2 


0.070 + 0.039 



Note. — This table reports the bolometric luminosities, the Gauss-Hermite expansion coefficients for the BF profile, the inferred 
inclination, and an estimate of the central black hole mass and accretion rate. 
Objects marked with a * have a polarization measurement in the NVSS catalogue. 



case of the BLR. 

A particularly important problem, involving the de- 
termination of AGN physical properties from emission 
lines, resides in the line profile asymmetries. Several fac- 
tors, such as partial obscuration, geometrical structure, 
or large scale non-virialized motions can produce asym- 
metric line profiles. Moreover, relativistic effects within 
the gravitational field of the SMBH give raise to asymme- 
tries, especially in the high velocity wings of the profile, 
which are included in t he calculations of the model by 
IChen fc Halpernl (119891 ). In order to assess how much 
the asymmetric component affects our estimates of the 
velocity field, we introduced an asymmetry parameter: 

K3 = hMBWEMiip), (20) 

expressing the relative contribution of the asymmetric 
component, with respect to the Gaussian component, in 
the profile of H/3 at its half-maximum level. As we show 
in Fig. 8, the asymmetric component gives a contribu- 
tion to the FWHM which rarely exceeds the 10% level. 
The most extreme cases, where the asymmetric compo- 
nent becomes larger than 20%, occur only in the range 
of very broad line emitting sources. Although this does 



not appear to be a general property of broad line ob- 
jects, it echoes the observation of larger asymmetries in 
objects where FWHMh/3 > 4000 km s~^, which is amo ng 
the features identified bv lSulentic et al. I (|2000l . l2006l ) in 
their distinction between Population A and B sources. 
Objects with the largest asymmetries are more problem- 
atic for the comparison with the reference model, used 
to calculate the equivalent velocity field in the BLR. 
However, comparing their masses and accretion rates 
with those of the other sources, we do not find system- 
atic differences that could suggest the need for specific 
model corrections in the asymmetric line emitters, as it 
is shown in Fig. 9. Here we see that the result of in- 
troducing the BLR inclination in our calculations is to 
remove the strong dependence of the accretion rate on 
FWHMh/3, that was commonly found with isotropic mass 
estimates. Instead, we are left with a much more com- 
plex situation, where, though a slight trend to measure 
lower accretion rates in broad line emitting sources is still 
present, it is not a universal condition. An inverse power 
law fit to the data yields Lboi/LEdd oc FWHMh/j"" "^^, 
considerably weaker than the old isotropic prediction of 
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Fig. 8. — H/3 asymmetry parameter distribution with respect to 
FWHMh/3. Filled symbols represent objects where the asymmetric 
component exceeds 10% of the Gaussian contribution. We plot as 
triangles objects that are affected by negative asymmetry, yielding 
blue shifted peaks and red shifted wings, and squares for positive 
asymmetry sources, having a red shifted peak with blue shifted 
wings. Large asymmetries characterize objects where fits with the 
reference model are more likely to be problematic. The bar in 
the upper left region of the diagram is a median estimate of the 
measurement errors. 

Lboi/LEdd FWHMh^"^. In particular, we do not 
observe dramatic excesses in the accretion rate of our 
sources, whose estimates are far below the correspond- 
ing Eddington limits. 

Most of the results achieved in this work depend crit- 
ically on the choice of our reference model, which leads 
us to conclude that the BLR has a flattened compo- 
nent, that is commonly seen at i < 20°. In the case 
of radio-loud sources, nearly face-on disk structures are 
likely to produce a radio jet oriented along our line of 
sight towards the object, and the resulting signal should 
be highly variable and polarized. Although we were not 
able to find any information about variability, some of 
the radio loud sources in our sample have been detected 
in the NRAO VLA Sky Survey (NVSS), which provides 
measurments of the radio flux and polariz ation at the fre- 
quency of 1.4 GHz (jCondon et al. 1119981 ).'^ We identify 
these objects in Table 2 and we compare the degree of lin- 
ear polarization with our inclination estimates in Fig. 10. 
Although the uncertainties are quite large, a significant 
degree of linear polarization is detected in many objects 
and it appears to be an averagely decreasing function of 
i. 

A comparison of our mass determinations with the old 
isotropic assumption allows us to study the properties 
of the geometrical factor within our sample. The sit- 
uation depicted in Fig. 11 clearly indicates that signif- 
icant effects, up to a factor ~ 30, should be expected 
and that they are more commonly observed in the range 
of sources with FWHMh/3 < 3000 - 4000kms-i. We 

^ Polarization data are available at 

|http://www.cv.nrao.edu/nvss/NVSSlist.shtml| 
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Fig. 9. — SMBH accretion rates as a function of FWHMh/3, with 
the same symbols as in Fig. 8. The cross in the upper right region 
of the diagram is the median uncertainty of measurements. We 
do not detect either dramatic accretion rate excesses or systematic 
trends associated to the line width, although higher accretion rates 
appear to be more common in the domain of narrow emission line 
sources. This conclusion seems not to be affected by asymmetric 
line profiles. 

find that the average value of the geometrical factor for 
black hole mass determinations based on FWHMh/3 is 
/ = 10.58 l b 7.70, margina.lly co nsistent with the result 
achieved bv lOnken et al. I (|2004f ). who gave / = 5.5± 1.9 
using the emission line dispersions. 

5. CONCLUSIONS 

In this work we investigated the shape of the emission 
line broadening function in the BLR of active galactic nu- 
clei. We used a technique based on cross-correlation and 
Gauss-Hermite line profile fitting, applied to the Balmer 
series, to infer the broadening functions and we compared 
them with the predictions of a structural model for the 
BLR. According to our results, we come to the following 
conclusions: 

• the line profile broadening functions carry much 
detailed information about kinematics of the BLR, 
which can be better understood by means of tech- 
niques exploiting the whole profile, rather than re- 
stricting on specific parameters; 

• the observed distribution of line profile kurtosis is 
consistent with the presence of a flattened compo- 
nent in the BLR, with a typical inclination i < 20°, 
though the actual values of i may depend on the 
adopted model; 

• some of the objects included in the sample have 
quoted measurments of linear polarization at radio 
frequencies, which averagely increase as the esti- 
mated BLR inclination approaches face-on; 

• correcting the SMBH mass and accretion rate es- 
timates for geometry reduces the anti-correlation 
among FWHMh/3 and accretion rate to a much 
weaker trend; 
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Fig. 10. — Degree of linear polarization at the radio frequency of 
1.4 GHz as a function of the inferred BLR inclination. The cross 
in the upper left corner gives the median uncertainty estimate. 
Polarization data are from the NVSS catalogue. 

• there are no particular indications, in our results, 
for a strong influence of the line profile asymmetries 
on the determination of the SMBH properties. 

Although this analysis may represent an advance in the 
problem of determining the role of the BLR geometrical 
factor, more questions should be answered, concerning 
how the two components combine in the observed line 
profiles. A precious contribution in this effort would 
probably result from the comparison of this technique 
with some independent way to estimate the black hole ac- 
cretion rate. Recent works suggested that this test could 
be possible, for example, with X-ray observations. The 
ability to constrain the BLR geometrical factor, then, 
could be applied to study the properties of black hole - 
host galaxy scaling relations with improved accuracy. 
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